Symmetry relations for spin-resolved exchange correlation kernels 
in soft magnetic layered systems 
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We first exploit the physical condition satisfying the symmetry relation of the "exact" spin- 
resolved exchange correlation kernel based on the "mixed scheme" in soft magnetic layered systems. 
The conditions are derived and examined by taking into account the field gradients of the magnetic 
moment as well as that of the electric moment. We also exploit the physical condition by means of 
deviation distribution function suitable for complex electronic structures with arbitrary £. 
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Density functional theory (DFT) has been a chief sup- 
port of the condensed matter theory for its practical pre- 
diction of various properties of many-electron system. It 
has resolved many properties such as the ground state 
energies, forces on the atoms or molecules, and their 
equilibrium positions Q. Its local density approxima- 
tions including density gradient corrections have gener- 
ally been used for their simplicity and accuracy 0, 
Unfortunately, there have been some exceptions for their 
accuracy. On top of that, semilocal approximations fail 
to describe the long range interaction between electrons 
whose densities do not or weakly overlap. As an alterna- 
tive, the "mixed scheme" has been proposed 0,0>ISEI0- 
In this scheme, the Coulomb electron-electron interaction 
is decomposed into the short range (SR) part and the 
complementary long range (LR) part. The key concept 
is then to use a semilocal approximation for SR part and 
a more appropriate many-body methods for LR part. 

In view of practical interests for spin transport prop- 
erties of soft magnetic layered systems, the resolution of 
topics of the weak dispersion or van der Waals forces is 
clearly important. Additionally, there has been grow- 
ing great deal of interests for the magnetic information 
stora ge syste m basically containing soft ferromagnets 
|sL ioL llfJ. Illj. These storage applications are based on 
the idea of the "spin-torque switching" in fcrromagnetic- 
normal-ferromagnetic metal hybrid structures. The idea 
has been proposed by Slonczewski and Berger 

as an attractive alternative of controlling the magnetiza- 
tion in place of the traditional way of modulating mag- 
netic fields. Considerable evidence for the concept has 
been accumulated by demonstrating that the magnetiza- 
tion can be switched or precessed in a soft ferromagnet, 
because a current exerts "spin-transfer torque" on a thin 
free magnetic layer, which is polarized in a thick "fixed" 
magnetic layer |1 0, El El El El Ei Such interests in- 
spired us to propose symmetry relations for the exchange 
correlation kernel (XCK) in soft magnetic systems involv- 
ing intricate boundary conditions and corresponding spin 
scatterings. 

As one of the central signatures in DFT, XCK can 



be useful for giving direct insight into the dynamics of 
spin interactions in multilayer system. XCK is a basic 
concept in describing many body correlation effects in 
an inhomogcneous electron liquid and is a vital element 
of the spin susceptibility. Spin-resolved XCK /™/(C) is 
defined by 



_ d 2 E xc (r,r';C) 
~ dn s <(r)dn s (r') 



(1) 



where E* c is the total XC energy functional in many 
electron system with spin polarization £. Here, n s (n s i) 
denotes the electron density in spin state s(s'). /* S C /(C) 
satisfies the symmetry relation 



mo = m-o 



(2) 



which is a key property of XCK and equally a fundamen- 
tal idea of importance in condensed matter. The specific 
symmetry relation 



/£(0 = /£ C (C) 



(3) 



also plays a significant role to investigate the situations 
in artificial composite structures However, prior to 
our previous work, there was no attempt to interpret 
XCK in terms of directly measurable quantities such as 
spin current I S ( S )- Although considerable studies have 
been devoted to investigate I S ( S ), theoretical approaches 
have substantially been limited to oversimplified multi- 
layer systems unsuitable for more realistic applications. 

The introduced key idea in this letter is the accu- 
rate relations guaranteeing the specific symmetry rela- 
tion in soft magnetic layered systems. The accurate re- 
lations are derived explicitly by including not only the 
field gradient of the electric moment but also that of the 
magnetic moment. We also explicitly show the relation 
which guarantees the specific symmetry relation by us- 
ing the solution of Boltzmann equation without neglect- 
ing anisotropic terms in systems of various dimensions. 
The condition satisfying the specific symmetry relation 
could be related to the spin density variation Vji s (ij and 



2 



then Vn s ( s -) to the variation of electrochemical poten- 
tial V/i s (j). Hence a general inhomogeneous multilayer 
system with arbitrary magnetization alignments can be 
explored easily by examining properties of XCK on the 
basis of 7 s(5 ) cx 

Density functional based on the "mixed scheme" is 
given by the sum of the kinetic energy T of a nonin- 
teracting particle system, SR and LR potential energies 
U(= Us + Ujj), and unknown SR and LR functionals 
E xc (= T xc + E£ c ); E[n{r)] = T[n{r)] + U[n{r)} + E* c in 
terms of one particle density in N-electron system 

n(r)=N ^2 J\^(rs,r 2 s 2 --- ,r N s N )\ 2 dr 2 ,--- ,dr N . 



Once SR (LR) part is given by T[n(r)]+[/ s(L) [n(r)]+£ xc , 
the complementary LR (SR) part is the difference be- 
tween the universal functional E[n(r)] and the SR (LR) 
part. In pair density theory giving a more improved 
ground state energy than in one particle density []]], 
the spin-summed pair density is given by n(r 'i,t*2) = 
N(N-i) J2 siS2 7 SlS2 (ri,r 2 ) where -f SlS2 (r 1 ,r 2 ) is the 
spin-resolved diagonal of the two-body reduced density 
matrix Q. The exact energy density functional is given 
by E[n(x)] — T[n(x)] + U with x — (ri,r 2 ). Here, 
U(= Us + Uh) consists of the external potential of a given 
pair and the interaction potential between particles form- 
ing two pairs at x*i — {rn,ri 2 ) and xj = {rji,rj 2 ) 0. 
Hence, the total XC energy functional T xc in an inter- 
acting system can be described in terms of pair density 
n(x) 

E* c [n(r)\ = T[n{x)\ - T[n{r)\ + U[n(x)] - U[n[r)}. 

E* c contains the unknown remaining kinetic term AT = 
T[n(x)] — T[n(r)] as well as the interparticle interaction 
potential. Here, AT denotes the difference between the 
exact kinetic term T[n(x)] in an interacting scheme and 
noninteracting counter part T[n(r)]. f££, the second 
derivatives of T xc with respect to spin densities n s (r) 
and rij(r') at a given pair position Xi — (r,r'), is now 
written in the form 



V r V r ,[AT[n(x)} + U[n(x)} - U[n(r)} 
Vn 5 (r)Vra a (r') 



(4) 



For the defined XCK by Eq.QJ, the symmetry relation 
fss(C) — fss(~~() 1S trivial in an inhomogeneous system 
of broken spin symmetry since majority and minority 
spins interchange their orientations and positions with 
the reversed polarization — £. That is, based on the as- 
sumption that the physical condition 



Vng(r') _ Vng(r) 
Vn s (r') Vn s (r) 



(5) 



is fulfilled (i.e., the ratios of spin density gradients are 
the same at two different positions r and r'), the specific 



symmetry relation 

9 2 T xc (r,r';C) 9 2 T xc (r, r'\ Q 



dn-g(r)dn s (r') dn s (r)dn s (r') 



f£(0 (6) 



is obtained trivially. So to speak, we explicitly show the 
specific symmetry relation by using XCK defined on the 
"mixed scheme". Hence it is elucidated that Eqs.© and 
© can be used to investigate soft magnetic systems with- 
out difficulty. Equation© thus can be useful to examine 
the adequacy of the specific symmetry relation in partic- 
ular for the present weakly interacting systems. 

The continuity equation for the probability density 
PM{t)[= PB(n s — Tls)] is written in form [lfij | 



dp M (t) 
dt 



-V-j M (t) 



(7) 



where the probability current density jijw (^) [= Mb(j's — 
js)] on Af-sphere is given by [j| 

3M{t) - pM{t)Mdet ~ DV PM (t) (8) 

where D is the diffusion constant in the case that diffu- 
sion constants of up-spin and down-spin are the same. 
Here, Mdct is the deterministic part of the Landau- 
Lifshitz (LL) equation for the evolution of a uniform 
magnetization M(t). The deterministic part is the sum 
of the conservative precession, dissipative LL damping, 
Slonczewski current-induced term and a neglected 
term in the current-induced part 0; 



Mdet = -7M x Hcons - yaMm x (m x 
—^fJMrh x (m x rh p ) 



^rn x (9) 
8?rM dt y ' 



where 7, M, a, and J are, respectively, the gyromag- 
netic ratio, saturation magnetization, LL damping con- 
stant, and an empirical constant with units of mag- 
netic field which is proportional to the current. Here, 
to = M/M and rn p is the unit vector along the mag- 
netization and the easy axis. In the first and second 
terms, H cons is the field around which M precesses as 
conservative and is related with the magnetic energy Em 
via p H cons = -^f-- The magnetic energy E M is the 
sum of the anisotropy energy and the coupling term with 
an external field |9|. When the last term (containing 
the mixing conductance g ss ) is neglected, the current- 
induced part in Eq.© is consistent with the result of 
Slonczewski Q. 

The current for up (down)-spin carriers can be written 

by 



3s(s 



M dc 



»(*)■ 



D s{s) Vn s{s) (10) 



where v is the mobility and D s fs) is the spin- up (down) 
diffusion constant. Eauation HlU|) is based on the fact 
that the drift velocity is proportional to the force on 
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the conduction electron. The force originates from the 
field gradients of its electric and magnetic moment (i.e., 

Here, the effective 

field H e ff includes the external field and the molecu- 
lar field coupled to the conduction electron via an ex- 
change constant. Hence, the spin-up (down) velocity 
v s(s) — —vE ± Mdct/2 is applied in Eq. l(TT)|) . because the 
difference between the velocities of up-spin and down- 
spin is given by v s — v§ = 2/jbV 'H c s/ k = M^ct from 
Eq.© 

The probability density pm often depends only on en- 
ergy, for instance, in a thermal equilibrium or a nonequi- 
librium steady state without LL damping a and current 
J. For this reason, the density pM(t) can be replaced 
with /9m (e, t) from the assumption that pm(s, t) depends 
only on energy. Still energy dependence can be different 
in distinct positions of sphere. Together with pjv4-(e,£), 
the spin-up (down) probability current density can be 
obtained by 



and I HE ^ 0, the condition which guarantees the specific 
symmetry relation is 



j s (s)0M) = f \j s (s) x dM] ■ m 



' D S(S, q 



( £ ,f){,/ EE T |[(«M. 
dn s{s) {e,t) HE 



^)/ HE -Jm p -/ M ]} 
(11) 



Here, energy integrals for the electric field J EE and the 
magnetic field 7 HE , and magnetization integral / are 
defined, respectively, by 

7 EE = <j)dMx E-rh 

7 HE = & dM x ff r , ona • m = I H c „n S dM 



(12) 



I M = IdMxM-m. 



On the assumption that the system is in steady state, 
the probability current density jM{t) vanishes. To get 
some msig ht in Eq.jni), we consider possible systems in 
steady state (i.e., Jm (*) = 0) an d investigate the specific 
symmetry relation in each feasible case. A steady state 
is reached if the condition j s {e) = js{s) is satisfied. After 
some straight forward algebra, we can write the relation 
corresponding to above condition j s (e) = Js(^) by 



Vn, 



Vn. 



(14) 



In this situation, the specific symmetry relation is triv- 
ially satisfied in normal region due to D s = D s . Another 
possible meaningful situation satisfying j s (e) — js(s) is 
that n s + rig = resulted in Vn s + Vji s = 0. As long 
as vI EE = and 7 HE = 0, the condition assuring the 
specific symmetry relation is given by 



Vn, 



-1. 



(15) 



In this case, Eq. l|15|) is in agreement with the previous 
condition [15J for the specific symmetry relation in non- 
metallic region of a homogeneous system with no space 
charge. Even though the term containing the mixing con- 
ductance_o ss in Eq.© is neglected as the case of Slon- 
czewski j8|, all the situations examined here satisfy the 
same conditions, Eas. (|14|) and (|15fl . assuring the specific 
symmetry relation in the similar way with no terms of 
g ss . As ever, these results are very useful to investigate 
spin related phenomena particularly in the normal metal 
region between two ferromagnets due to the condition 
D s = D s . 

Regarding realistic systems such as quasi-two or three 
dimensional multilayer and complicated electronic struc- 
tures, corresponding spin-dependent scattering need to 
be considered. Consequently, the neg lected anisotropic 
terms and coordinate variables fl9ll20l| have to be consid- 
ered in electron distribution functions. For this purpose, 
we examine the spin current density given by |17| 



Js(s 



d 3 vvgf s \s,r) 



(16) 



in terms of deviation distribution function (e, r) cor- 
responding to the energy e and velocity v. A possible 
solution is given by 



9 s v (S) (e,r) 



eE 



dfoje ) 
de 



(17) 



Mo I 



nE^ Mn^Va j =vl^(n s {e)-n- s {e)) 



-2[(aM - 2^)1™ - Jm p ■ I M ](n s (e) + n s (e))(13) 



Equation (|13J) leads us to following feasible situations 
in the system. First consider the case that up-spin and 
down-spin densities are the same. As long as the ratio 
of the work of the current induced torque to that of the 
LL damping rh p ■ / M // HE is equal to (aM - 2§§P) /J 



for a system with all diffusive boundary effects driven by 
an applied electric field. It is valid in systems of various 
dimensions without neglecting anisotropic terms. Here, 
the effective mean free path is given by ^eff r ) = 

1 — cxp(^jpy^y) corresponding to an electron 

with normalized velocity v passing through a position r 
in the direction t*o — r with ro on the cylindrical bound- 
ary. The intrinsically anisotropic mean free path A s ( s )($) 
is given by Aq^^I — a s< - s 'cos 2 i9 — 6 s ( s 'cos 4, i9] where $ is 
the angle between the magnetization M and velo city v. 
When we compare Ea. (|17f) with another solution 18], 



v ■ eE\ s c ^ ] (■&, r)=p- fi s(S) (r) + v- eEX 



s(s). 



s(S) 



(18) 
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where /2 is the equilibrium electrochemical potential. Ac- 
cordingly, the relation convincing the specific symmetry 
relation is given by 



N s (e F ) V{v -eE[l- A* ff (tf, r)] } - eE 



N- s (e F ) V{t) -eE[l- Af ff (tf,r] } - eE 



Vn, 



(19) 



Eciuation (|19|) is based on the relation given, in the pres- 
ence of an electric field (E — -V3>h bv |20| 

Vn s ( s -) = eN s ( S) (e F )[V Hs {g ) - eE}. 

Here N s ^ S )(e F ) is the spin-up (down) density of states 
at the Fermi level in a highly degenerate system. In or- 
der to generalize Ea.Q19|). we compare <7« (e, r) given by 
Ea. (|18|l with another solution |2l| obtained by using the 
path integral along v. The gradient of the electrochemi- 
cal potential is then obtained by 

V/ia(s)(r) = -y[S s{s)l3 (r,r')g^(v,r')Sl s{s) (r 7 r')] 

-V / dl"S s(s)(3 (r,r")v ■ eE(r")Sl s{s (r,r")^ 

(20) 



+V(v-eEX s {s) 



where the spinor propagation factor is given by S(r, r') — 
P r '^ r exp(—^J r , r ,sdl"^(r")) with the inverse mean- 
free path operator £(r"). Here, P r '~, r is the path order- 
ing operator along T(r, r 1 ) which stands for the oriented 
straight path that starts at point r' and ends up at r. 
If we take the limiting case of sufficiently large value of 
\r — r'\, Eq. ((2U|) is written by 

V/x s(s -)(r) = -V/ r(r)0 dl"[S sm (r,r")v ■ eE(r") 

Sl^ry^+Viv-eEX^). (21) 

The first term on the right hand side in Ea. f21|) 
corresponds to — V(« • eEX^g (v, r)) in Ea. lfT^I as- 
suring the specific symmetry relation. Hence the 
condition assuring the symmetry relation is written by 

N s (e F ) Vli-eE\°-f r(r rl) dl''S s0 (ry')v-eE(r'')S\ s (ry')^]-eE 



Vn, I 
V)is I r 

Hence, 



NileF) Vlv-eE\^J r(T , tT . l) dl"S S 0(r,r")v-eE(r")S!, s {r,r")^-}-eE 

= , in agreement with Ea. lfT^|l . 

Ea . l|19|l and above related condition can be 
used to account adequately for experimental results. 
Once some quantities such as v, E and the angle d at 
distinct positions are known, we can check the proper- 
ties of XCK with ease. Besides, we can determine the 
relevance of theoretical and experimental consequences, 
particularly in the region satisfying the specific symme- 
try relation such as a homogeneous nonmagnetic system. 
For example, in case of choosing correct factors except 
an inappropriate X^g \ Eq. l|19|) is not satisfied in trivial 



situations and proper modification for components such 
as a s ^ and fe 5 '^ of A*g would be done to be in agree- 
ment with the relation. Also experimental values can be 
investigated as well and be compared with the modified 
calculations. 

Here, we have shown the conditions confirming the spe- 
cific symmetry relation in current soft magnetic layered 
systems. First, the specific symmetry relation has been 
elaborated by defining XCK explicitly on the "mixed 
scheme" compared with our pervious work [l5| . Then 
the conditions are derived by taking into account the 
field gradients of the magnetic moment in addition to 
that of the electric moment. At the same time, spe- 
cial attention is paid to the fact that when we consider 
the neglected anisotropic term in Boltzmann equation, 
it gives deviation from generally accepted definition of 
the electrochemical potential. For the first time, we have 
thus given the physical condition by dealing with devia- 
tion distribution function suitable for complex electronic 
structures in addition to one dimensional systems. 

In conclusion, we have proposed a new method to con- 
firm the relevance of studies of the spin-related phenom- 
ena in soft magnetic layered systems of various dimen- 
sions by taking the physical condition = 
into account. Once some quantities in the condition are 
given, the properties of XCK can easily be checked the- 
oretically and experimentally particularly in the region 
where the situation for convincing the specific symmetry 
relation is known. Furthermore, excellent agreement of 
theoretical and experimental results can be induced by 
modifying factors or experimental conditions in accord 
with the physical condition. 
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